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Abstract 

The effect of the Coulomb mteraction on the phase diagram of finite nuclei is studied within 
the Canonical Thermodynamic Model. If Coulomb effects are artificially switched off, this model 
shows a phenomenology consistent with the liquid-gas phase transition. The inclusion of Coulomb 
does not significantly affect the phase diagram but it drastically modifies the nature and order 
parameter of the transition. A clear understanding of the phenomenon can be achieved looking at 
the distribution of the largest fragment produced in each fragmentation event. Possible connections 
with experimental observations are outlined. 
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I. INTRODUCTION 



Nuclear multifragmentation in intermediate energy heavy ion collisions continues to be 
a topic of intense interest. In particular it has been recently proposed [l| that a bimodal 
behavior in the distribution of the heaviest fragment produced in each fragmentation event 
can be an experimentally measurable signature of a first-order transition. Different data- 

nnn 

sets have been reported to confirm such an expectation [2|, |3|, |^ and a lively debate on its 
ultimate physical interpretation was raised . To progress on these issues, here we 

report on calculations using the canonical thermodynamic model (CTM) [9]. A study of the 
thermodynamic properties of the model will show that indeed the bimodality signal can be 
associated to the finite system counterpart of a first-order transition in the framework of 
our model, but that the nature of such a transition is deeply modified by the presence of 
Coulomb with respect to ordinary liquid-gas. 

Our ultimate aim is to confront the model with some data which exist from quasi- 
projectile (QP) fragmentation of Au on Au 0, [ll| but we will also do, in the beginning, 
some calculations for hypothetical nuclei with Coulomb interaction switched off. The rea- 
son for this preliminary study is that the role of the Coulomb interaction on the observed 
phenomenology can be clearly spotted by such a study, interaction. 



II. BASIC EQUATIONS OF CTM 

r 

The mathematical machinery for calculating with CTM has been described before [9] but 
to establish our notation and also list the parameters used in this calculation we need to go 
over some details. The dissociating nucleus breaks up into clusters. Each cluster is specified 
by two indices a=number of nucleons and j=number of protons. The canonical partition 
function for the fragmenting source with A nucleons and Z protons (neutron number N = 
A — Z) dX a, given temperature T is given in our model by 

Qa,z = Y.X{^ (1) 

Here the sum is over all possible channels of break-up which satisfy the conservation laws; 
ujaj is the partition function of one composite with nucleon number a and proton number j 
respectively, and .,■ is the number of this composite in the given channel. The one-body 
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partition function cuaj is a product of two parts: one arising from the translational motion 
of the composite and another from the intrinsic partition function of the composite: 

Uaj = ^(27rmaT)3/2 x Za,j{int) (2) 

Here ma is the mass of the composite and Vf is the volume available for translational motion; 
Vf will be less than V, the volume to which the system has expanded at break up. We use 
Vf — V — Vq , where Vq is the normal volume of A nucleons and Z protons. We will shortly 
discuss the choice of Zaj{int). The freeze-out density in unit of normal nuclear density is 
p/po = Vo/V. 

The probability of a given channel P{na,j) = -P(t^i,Oj "^i,!; "^2,1 n,aj ) is given by 

Pina,) = 7^U^ (3) 

The average number of composites with a nucleons and j protons is seen easily from the 
above equation to be 

/ \ QA—a,Z—j /,\ 

(riaj) ^ u;a,j—pz (4) 

The constraints A = x Uaj and Z = J2j >^ f^a.j can be used to obtain different looking 
but equivalent recursion relations for partition functions. For example 

Qa,Z = ^Y^3^a,jQA-a,Z-j (5) 

We now give the choice of Zaj{int) used in this work. The proton and the neutron are 
fundamental building blocks thus zifi{int) = zi^i{int) = 2 where 2 takes care of the spin 
degeneracy. For deuteron, triton, ^He and ^He we use Zaj{int) — {2sa,j + 1) exp{—Peaj{gr)) 
where /3 = 1/T,eaj{gr) is the ground state energy of the composite and {2sa,j + 1) is the 
experimental spin degeneracy of the ground state. Excited states for these very low mass 
nuclei are not included. For mass number a — 5 and greater we use the liquid-drop formula. 
For nuclei in isolation, this reads 

Za,j{int) = exp -[Woa - a{T)a^/' - - ^ + 6 

The expression includes the volume energy, the temperature dependent surface energy, the 
Coulomb energy, the symmetry energy and contribution from excited states in the continuum 
since the composites are at a non-zero temperature. 



In using the thermodynamic model one needs to specify which composites are allowed in 
the channels. For mass number a=5 we include proton numbers 2 and 3. For mass number 
a=6 we include proton numbers z=2, 3 and 4. For higher masses we have followed this 
procedure. The liquid-drop formula allows one to define drip-lines for a given a. For a > 7 
we include all nuclei within drip-lines. This choice allows us to use the same criterion in all 
cases studied, i.e., when the Coulomb is switched off or half turned on or fully turned on 
(required to compare with actual data). We can not prove if our choice of composites is the 
best one to use from the point of view of principles. But at least it is well defined. Some 
study was made in {l^ on the effects of changing the width of the ridge of the nuclei used 
in computation of properties we seek. Here we have stuck to one prescription. 

The Coulomb interaction is long-range. The Coulomb interaction between different com- 
posites can be included in an approximation called the Wigner-Seitz approximation. We 
incorporate this following the scheme set up in jo, [l3|. This requires adding in the ar- 
gument of the exponential of Eq.(l7l) below a term K,-^^{Vo/Vy^^ . Defining Rj = ^ the 
average energy of the system is given by (E) = + where for a > 4 we have 

e,^^ = lT + a{-W, + T^/e,)+a{T)a^ri + K^,[l.Q-{Vo/VYl^] + s^-^ 
For a < 4 we use Caj = |T -|- Cajigr) — K-^{Vo/Vy^^ . We label as (E*) the excitation 
energy: (E*) = (E) — E{gr) where E{gr) is calculated for mass number A and charge Z 
using the liquid-drop formula. The pressure in the model can be shown to be simply 



III. THE LARGEST AND THE SECOND LARGEST CHARGE IN EVENTS 

As recalled in the introduction, the size of the largest fragment (or equivalently its atomic 
number Zi) is an especially interesting observable in the multi-fragmentation problem. Not 



only it is an experimentally accessible quantity in exclusive experiments 



14| , but it is known 



to provide an order parameter of fragmentation transitions for a large class of equilibrium 



as well as out of equilibrium models [15|]. The second largest fragment ^2, though not so 
important from the theoretical viewpoint, has been also extensively used to characterize the 
topology of fragmentation [igI, 17|. 



in an event, by Z2. Calculations of Z2 in CTM have not been done before and in order to 
derive the expression, it is advantageous to derive the one for zi first. 
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There is an enormous number of channels in Eq.(l). Different channels will have 
different values of zi. For example there is a term -^-|f- this channel the 

highest value of Z is 1 and thus zi is 1. The probability of this channel occurring 
is (from Eq.(l3])) Q^^r^r- The full partition function can be written as Qa,z = 
<5A,z(i^i,05 "-^1,15 "^2,1, ^3,1, i^3,2 i^aj )• If WG coustruct a Qa,z where we set all tu's ex- 

^Z 

cept cui^o and Ui^i to be zero then this Qa,z{^i,o, t^i.i, 0, 0, ) = -^^r and this has 

Zi = 1. Clearly, Qa,z{^^i,o, 0, 0, ) is contained in and is a small part of the full 

partition function Qa,z = <5A,z(t^i,o, ^^2,i, tus.i, ^^3,2 ^^aj )• 

It will be convenient to introduce a shorthand notation. Except when confusion may 
arise, we will write Uj to collectively mean all of uj{a,j) where j is fixed but the sum over a 
runs over the allowed range. As noted above, the full partition function denoted by Qa,z is 

Qa,z = Qa,z{^i,q, ^^1,1,^^2,1, 1^3,15 ••••) 

= Qa,z{^Q, ^1,^2, ^z) (7) 

We are now ready to write down a general formula. Let us ask the question: what is 
the probability that a given value Zm occurs as the maximum charge? To obtain this we 
construct a Qz,n where we set all values of tu^ = when z > Zm- Call this QA,z{zm)- 
Then QA,z{zm)/QA,z (where Qz,n is the full partition function with all the tu's) is the 
probability that the maximum charge is any value between 1 and z^- Similarly we construct 
a QA,z{zm — 1) where ujz is set at zero whenever z > Zm — ^- The probability pi that zi is 
Zm is given by 

Mz,) ^ Mz,^ = z,) = Q"'"' - Q^'^'" - ^ (8) 
The average value of zi is 

{z^) = Y.^ipM) (9) 

and the rms of normalised zi is 

rms = ^Y.pM){{^iIZ?-{{z,)IZY) (10) 

We now turn to the calculation of 2:2- A given charge Z2 can be the second highest charge 
when (a) there is at least one particle in Z2 but also just one particle in a charge state zi > z^ 
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or (b) there is no no particle with charge z > Z2 but Z2 occurs at least twice. The partition 
function Qa for the case (a) is 

Qa = ^ [QA-a,Z--zi{z2) " QA~a,Z-ziiz2 - 1)] (11) 

For a fixed zi the sum over a is over all nuclei within drip-lines. The sum then goes over all 
Zi's larger than Z2- 

For the case (b) we have 

Qb = Qa,z{^2) - Qa,z{z2 - 1) - Xl^«.^2<5A-a,Z-22(^2 " 1) (12) 

a 

Finally ^2(^2)5 the probability that the second largest has the value Z2 is given by 

P2(^2) = [Qa + Qb]/QA,Z (13) 

The average value of Z2 results: 

(^2) = E ^2^2(2:2) (14) 

IV. COULOMB EFFECTS ON THE FRAGMENTATION TRANSITION 

For a hypothetical nucleus of mass A = 150, atomic number Z = 75 and with no Coulomb 
force, we plot the functional relation between pressure and total density in the upper left 
part of FiglH A few comments about EoS in the model are in order. 

The choice of a symmetric system was done to avoid mixing density and isospin effects, 
as we now explain. It is well known that phase transitions with two conserved charges (here: 
proton number Z and neutron number) can produce fractionation, i.e. different concen- 
trations in the different phases. Then to spot the phase diagram of such a system the two 



densities have to be varied independently 
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201 1 . A projection of the two-dimensional 



equation of state on a specific for example the P — p correlation at fixed N/Z, 

may be misleading, being continuous even in the occurrence of a discontinuous (first-order) 



transition [2l|]. This complication does not arise for symmetric matter, where fractionation 
disappears, the direction of phase separation follows the isoscalar density, and the whole 
information about the phase diagram is contained, as for one-fiuid systems, in the P — p 
equation of state. A first-order phase transition in this representation is unambiguously de- 
fined by the presence of a back-bending. Such behavior is very well known in the framework 



of the mean-field theory where it reflects the instability of homogeneous system with respect 
to phase separation. There, the appearence of a spinodal region reflects the inadequacy of 
the model. One needs to do a Maxwell or Gibbs construction before any correspondence 
with experiment can be established. It should be pointed out that the EoS in models like 
CTM are very different from those in the mean-field approximation, although a superficial 
look at Fig{T] would suggest this. There is a spinodal region in Fig{T]just as there is one 
in mean-field models, and in both cases they reflect phase separation, to see the effects of 
iso-spin in the model. The two EoS are quite similar. However, the backbending in CTM is 
an expected feature in the exact evaluation of the thermodynamics of a finite system which, 
at the thermodynamic limit, exhibits the discontinuity characteristic of a first-order phase 



transition 



22l | . The pressure decrease with density is physical, and it arises because of finite 



particle number imposed in the calculation. For earlier discussions on this see |23l . |2J]. In 
particular, in a simpler version of CTM it was explicitly demonstrated that the backbending 
disappears as the system grows and becomes a straight hne with zero slope as expected in 
a first-order transition 251] . The height of the straight line, that is the value of the transi- 
tion pressure for a given temperature, is given by the average pressure in the backbending 



region 



26|. 



Because of the preceeding discussion, we can conclude from Fig{T] that CTM presents 
a first-order transition between a dense, liquid-like phase and a dilute, gas-like one. The 
resulting phase diagram is reported in FiglH where the average pressure in the backbending 
region has been reported as a function of the temperature. The qualitative similarity with 
liquid-gas is apparent. It has been usual to limit the EoS in the in CTM to p/po < 0.5. The 
argument is that at higher densities the approximation of replacing the effect of the residual 
nuclear interaction between composites by just a constant excluded volume will begin to get 
worse. Here however both in FigsiT] and [2] the curves go beyond this. 

We consider now bimodality in the value of pi{zi) as a function of zi where pi(^i) is 
the probability that the charge of the largest cluster in an event is zi. The significance of 
a bimodal distribution in a finite system as a signature of first-order transition has been 



extensively discussed in the literature 126|, 



271]. Here we do the calculation at a given 



density p for different temperatures to locate the temperature where the two peaks of the 
bimodal distribution are nearly the same height. Having located the density and temperature 
we can then plot this as a point in the p — p plane. At the same density we can also 
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calculate the specific heat as a function of temperature and find the temperature where 
Cy maximises. This gives us another point in the p — p plane. Repeating the calculations 
for different densities we obtain two curves in the p — p plane. A maximum in the value 
of Cjj is a generic signature of a phase transition with finite latent heat rounded by finite 
size effects. In the context of nuclear physics, it was taken as a possible signature of the 
transition associated to multifragmentation a long time ago 28 1. 

The remarkable thing is that without the Coulomb interaction the line representing bi- 
modality and the maximum of coincide and fall in the backbending region. Both curves 
correspond to maxima of fluctuations: a peak in Cy = (t'e/T'^ means maximal energy fluc- 
tuations, and fluctuations in the size of the largest cluster are maximum at the bimodality 
point. In a statistical ensemble where the freeze-out volume would be free to fluctuate under 
the effect of an external pressure, it is easy to show that the system would also present a 
bimodality in the volume distribution and a maximum in volume fluctuation at the pressure 
corresponding to the backbending |26l. l30|. 

The coincidence of all these fluctuation signals agrees with the expectation that the 
fragmentation transition is very close to ordinary liquid-gas. Increasing the available volume 
at constant temperature, the system experiences a transition from configurations dominated 
by a single huge cluster at low excitation energy, to a highly fragmented pattern at high 
excitation energy. This transition can be classified as first-order in the sense that it would 
become discontinuous at the thermodynamic limit. In more technical terms, we can say 
that energy, particle density and zi are all order parameters of the observed transition. In a 
system as small as A = 150 of course all observables are continuous, however the transition 
point can be uniquely determined by the location of the jump (for p and zi) or the fluctuation 
peak (for E*) of the different order parameters. 

This is shown by the phase diagram of FiglJl The three different observables lead to a 
consistent estimation of the transition line. 

The situation changes with Coulomb force. We introduce a strength parameter Xc for the 
Coulomb force: Xc=0 means no Coulomb and Xc = 1 means the full Coulomb. 

Figs. [1] and [2] report the results on the isotherms, fluctuation loci and phase diagram as 
the strength of the Coulomb force changes from Xc = through intermediate values to its 
full value Xc = 1. 

Let us look at the P—p EoS first. We can see that the inclusion of the Coulomb interaction 
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results in a shift of the isotherms, with a gas-hke branch appearing at higher density respect 
to the uncharged case, and a lower transition temperature for a given pressure. Moderately 
charged systems up to Xc = 0.64 still show a phase diagram qualitatively similar to the 
uncharged case, but for the physical system Xc = 1 the isotherms monotonously grow and the 
phase transition has disappeared. These features are in qualitative agreement with previous 
works {29!. I30I . as well as with the well-known expectation that the transition temperature 



should vanish approaching the drip-lines 



3ll. 



The information given by the distribution of the largest fragment closely follows these 
thermodynamic findings. The bimodality line falls in the middle of the spinodal region for 
all values of Xc, and bimodality disappears together with the phase transition in the fully 
charged Indeed the system Z = 75, A = 150 is beyond the proton drip-line of 

the model, meaning that huge clusters close to the size of the source are unbound even at 
T=0. 

The situation is different for the curve of maximum c^- This curve can be defined even 
for Xc = 1, where the system is fragmented at any temperature and no phase transition is 
observed. This shows that a peak in c„ is a necessary but not sufficient condition for a first- 



order transition 



321 ]. When the transition is present, Figj2] shows that the observable leads 



to a definition of the phase diagram consistent with the isotherms and with the bimodality. 
However, we can see from Fig{T] that as the strength of the Coulomb force changes from 
Xc = through intermediate values to its full value, the splitting between the indicators of 
phase transitions progressively increases. Specifically, energy fiuctuations peak at a density 
lower than the one one corresponding to the spinodal region in the charged system. The 
fiuctuation energy peak is thus obtained within a pure phase. 

Of course experimentally we only see Xc=l. theoretically. As discussed above, the disap- 
pearence of the transition for Xc=l is due to the fact that the considered system is beyond 
the proton drip-line. The more neutron rich system A = 150, Z = 60 has the same size and 
the same Coulomb energy as our artificial A = 150, Z = 75 with Xc = 0.64. This nucleus is 
not isospin symmetric, meaning that the P — p correlation is not enough to trace the phase 
diagram of this system and a complete thermodynamic study would need the analysis of the 
whole two-dimensional P{pn,Pp) equation of state. 

However Fig|3] shows that isospin effects have a small influence in the P — p correlation 



of CTM. This is at variance with mean-field models for infinite nuclear matter 



20 



2l|. This 
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is not so surprising recalling that in CTM the pressure is simply proportional to the total 
fragment multiplicity j^, and this latter is only slightly affected by the isospin. Thus the 
P — p isotherms of the asymmetric system A = 150, Z = 60 are very close to the isotherms 
of the corresponding symmetric system shown in FiglH once Xc is chosen such as to have the 
same Coulomb energy as in the physical system. The same is true for the bimodality and 
Ct, signal, which are typical isoscalar observables strongly affected by Coulomb but not very 
sensitive to isospin. The preceeding discussion implies that the modification of the phase 
transition observed for < Xc < 1 in Figs IT][2] may be observed in a physically existing system 
like A = 150, Z = 60. This latter is close to the fragmention sources studied experimentally 
in refs.|2|,y,y,Q. 

The splitting of the two fluctuation signals observed in FigO suggests that the nature of 
the phase transition in charged systems may be different from that of uncharged ones. In 
particular in a transition with non-zero latent heat like liquid-gas the heat capacity should 
peak at the transition point, at variance with our results. This means that energy does 
not seem a good order parameter of the transition observed in CTM for charged nuclei. 
This last flnding is very close to the study of ref.{3o|. In that work it was shown that in 
the MMM fragmentation model in the isobar canonical ensemble the two bimodality peaks 
indicating the phase transition occur at very different Coulomb energies for heavily charged 
systems, but correspond to very similar excitation energies. If two phases are associated to 
the same energy, a mixing of the two does not induce any energy fluctuation, and the energy 
fluctuation at coexistence is the same as the energy fluctuation of the pure phases. Therefore 
we may expect that in the transition shown by CTM the two phases should correspond to 
similar excitation energies, or in other words that the latent heat of the transition should 
vanish with Coulomb. 

To progress in the understanding of the observed phenomenon, we show in FigH] the 
distribution of the largest zi and second largest Z2 fragment at the bimodality point for the 
same system as in Figj3l Not only the transition temperature is lowered by the inclusion 
of the Coulomb interaction as already observed, but the shape of the distribution is very 
different in the presence [xc = 1) or absence {xc = 0) of Coulomb. In the absence of Coulomb 
the two peaks are similar in shape and close to what is expected for the liquid-gas phase 
transition as depicted by the Lattice Gas Modelj^: a liquid-like solution with a dominant 
cluster exhausting about 75% of the total mass, while the gas-like solution appearing at 
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the same temperature consists of a much more fragmented configuration where the largest 
fragment is about 25% of the total size, and the second largest is of comparable size. The 
situation is completely different for the physical case where Coulomb is accounted for. In the 
presence of Coulomb the high-density solution corresponds essentially to the initial nucleus 
excited in its internal states. Such configurations give rise to what in nuclear physics is 
called an evaporation residue (recall secondary decay is not accounted for in CTM). The 
low-density solution is very different from the picture of a nuclear gas: the largest cluster 
is peaked at zx = Zjl and the second largest has a very broad distribution ranging from 
zi = Z/12 to Z2 = Z/2, and is close to the phenomenology expected for hot asymmetric 
fission. 



V. SUMMARY 

In this paper we have studied the effect of the Coulomb interaction on the fragmentation 
transition which is observed at finite temperature in the Canonical Thermodynamic Model. 
The typical behavior which is expected for a finite system counterpart of the liquid-gas 
phase transition is only observed if the Coulomb interaction is artificially switched off. The 
transition temperature decreases with increasing Coulomb energy as observed already in 



many other models 



a, 



30 



3l| . More interesting, the fragmentation pattern associated 
to the transition changes completely in the presence of Coulomb. In the spinodal region 
defined by the backbending of the p — p isotherms the distribution of the largest fragment 
is bimodal. The two dominant fragmentation patterns defined by the two peaks at the bi- 
modality point do not correspond to a liquid-gas phenomenology but are close to a transition 
from evaporation to asymmetric fission. 

CTM was recently shown to produce results which are close to another very successful 



model of nuclear fragmentation, the Copenhagen SMM[33|]. Because of that, we think that 
the presented results are not specific to our model but should be characteristic of any model 
of fragmentation in statistical equilibrium. 

The distribution of the largest fragment as a possible signature of a fragmentation tran- 
s,«on . expensively .ud.ed expe.,„e,.aUyE. flta i. quas.-p.oject.le .a„„e,.at,o. of 



Au+Au collisions. The comparison of CTM with experimental bimodality data will allow 
to progress on the interpretation of the transition observed in the data which is presently 
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the object of intense debate [5|, lO, l7|, l8|. This will be the subject of a forthcoming paper. 
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FIG. 1: EoS {p — p) for hypothetical nuclei with a varying strength of the Coulomb force Xc from 
zero (upper left) to the physical value Xc = I (lower right). The temperatures (in MeV) of the 
isothermals are indicated. Along the solid line the two peaks in pi{zi) as a function of zi are of 
the same height; pi{zi) is the probability that the maximum charge is zi. The dotted line denotes 
the line where the specific heat Cy is maximum. The two lines coincide for Xc = 0. 
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FIG. 2: Transition pressure as a function of the temperature for an hypothetical isospin symmetric 
A = 150, Z = 75 nuclear system calculated with the three different indicators of the transition: 
the average pressure in the backbending region (squares), the curve along which the two peaks in 
Pi{zi) have equal heights (circles), and the curve along which the c„ is maximum (crosses). The 
strength of the Coulomb force is varied from Xc=0 to Xc = 0.64. 



15 




FIG. 3: P — p correlation for a nuclear system with Z = 60,A = 150 without (xc = 0) and 
with {xc = 1) inclusion of the Coulomb force. The temperatures (in MeV) of the isothermals are 
indicated. Along the solid line the two peaks in pi{zi) as a function of zi arc of the same height; 
Pi{zi) is the probability that the maximum charge is zi. The dotted line denotes the line where 
the specific heat Cy is maximum. The two lines coincide for Xc = 0. 
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FIG. 4: Probability pi{zi) as a function of zi when the bimodal distribution has nearly same 
heights for the two peaks. The system is ^4=150 and Z = 60. The strength of the Coulomb force 
Xc has been progressively changed from (no Coulomb) to 0.25 to 0.5 to 1 (normal strength). Note 
that the shapes of the curves change very significantly. The calculation is done at a fixed density 
p/po=l/3 and the temperature varied till the bimodal distribution appears. The temperatures in 
MeV are indicated. 
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